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Abstract 

We present a simultaneous analysis, within an impact parameter dependent saturated 
dipole model, of exclusive diffractive vector meson (J/tp, 4> and p) production, deeply 
virtual Compton scattering and the total "y*p cross section data measured at HERA. 
Various cross sections measured as a function of the kinematic variables Q 2 , W and t are 
well described, with little sensitivity to the details of the vector meson wave functions. 
We determine the properties of the gluon density in the proton in both longitudinal and 
transverse dimensions, including the impact parameter dependent saturation scale. The 
overall success of the description indicates universality of the emerging gluon distribution 
and proton shape. 

1 Introduction 



Exclusive diffractive processes at HERA, such as exclusive vector meson production or deeply 
virtual Compton scattering (DVCS), are excellent probes of the proton shape in the perturbative 
regime. Several investigations have already shown that these processes can be well described 
within a QCD dipole approach with the vector meson wave functions determined by educated 
guesses and the photon wave function computed within QED; see, for example, Refs. [1-11]. 
It was also pointed out some time ago that the exclusive vector meson and DVCS processes 
provide severe constraints on the gluon density at low-x [12-21]. 

The vector meson and DVCS processes are measured at HERA [22-31] in the small- x 
regime where the behaviour of the inclusive deep-inelastic scattering (DIS) cross section, or the 
structure function F 2 , is driven by the gluon density. The dipole model allows these processes 
to be calculated, through the optical theorem, from the gluon density determined by a fit to the 
total inclusive DIS cross sections. Usually, it is assumed that the evolution of the gluon density 
is independent of the proton shape in the transverse plane. The investigation of Kowalski 



and Teaney (KT) [1] has shown that the Gaussian form of the proton shape, implied by the 
data, has implications on the emerging pattern of QCD evolution and saturation effects. The 
interplay of saturation and evolution effects was first investigated by Bartels, Golec-Biernat 
and Kowalski [32], where it was found that the total inclusive DIS cross sections, or F2, can 
be described either by strong saturation and weak evolution or by strong evolution and weak 
saturation effects. The investigation of Ref . [1] , which took into account also the proton shape 
in the transverse plane, concluded that saturation effects are substantial in the proton centre, 
but that the Gaussian form implies that a large contribution to the cross section has to come 
from the outskirts of the proton, where the gluon density is diluted. Hence the evolution effects 
have to be strong and play an important role. An alternative approach to determining the 
impact parameter dependent gluon distribution, based on a two-Pomeron model, is discussed 
in Refs. [33,34]. 

Another important result of dipole model investigations is that a wide variety of DIS data 
can be described with only a few assumptions. The investigations of Refs. [32,35-38] show that 
the inclusive DIS cross section can be described together with the inclusive diffractive DIS cross 
section. Moreover, in Ref. [1] it was shown that the inclusive DIS process can be described 
together with inclusive charm production and exclusive diffractive J/ip photoproduction. This 
description preserves also the main properties of the inclusive diffractive DIS cross section [39]. 

In this paper we will extend the analysis of Ref. [1] and show that the same minimal set of 
assumptions allows the description of a much wider set of recently measured data on exclusive 
J/ip, <p and p photo- and electroproduction and also the DVCS process. The cross sections for 
these processes have been measured as a function of the photon virtuality, Q 2 , the 7*p centre- 
of-mass energy, W, and the squared momentum transfer, t. In addition, for vector mesons the 
ratios of the cross sections for longitudinally and transversely polarised incoming photons have 
been determined as a function of Q 2 . 

To perform the analysis we use an impact parameter dependent saturated dipole model in 
which the gluon density is determined by a DGLAP fit to the total inclusive DIS cross sections. 
The wave function of the virtual photon is known from QED and the proton and vector meson 
wave functions are assumed to have a Gaussian shape. The parameters of these Gaussian 
distributions are easily determined from data. The results are compared to numerous data 
distributions provided by the HERA experiments. In this framework the W distributions are 
mainly sensitive to the square of the gluon density and the Q 2 distributions and ul/ut ratios 
to the properties of the vector meson wave functions. The proper choice of the wave functions 
is also confirmed by the agreement of the predicted size of the cross sections with data. In the 
dipole model the absolute normalisation of the vector meson cross sections follows from the 
optical theorem. 

The t-distributions determine the area size of the interaction region, B D . The parameter 
B D is obtained by making a fit to the t-distributions of the form da/dt oc exp(—B D \t\). For 
scattering of very small dipoles Bo is connected to the proton radius R p via Bo = Rp/S. 
However, for larger dipoles the size of the interaction area depends not only on the proton 
radius but also on the size of the produced vector meson or real photon, which we take into 
account following the work of Bartels, Golec-Biernat and Peters (BGBP) [40]. This allows the 
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Figure 1: The elastic scattering of a virtual photon on a proton in the dipole representation. 

data for all vector mesons and DVCS to be described using a unique Gaussian proton shape, 
independent of the produced final state. 

2 The dipole model 

In the dipole model, deep inelastic scattering is viewed as the interaction of a colour dipole, 
that is, mostly a quark-ant iquark pair, with the proton. The transverse size of the pair is 
denoted by r and a quark carries a fraction z of the photon's light-cone momentum. In the 
proton rest frame, the dipole lifetime is much longer than the lifetime of its interaction with the 
target proton. Therefore, the elastic j*p scattering is assumed to proceed in three stages: first 
the incoming virtual photon fluctuates into a quark-antiquark pair, then the qq pair scatters 
elastically on the proton, and finally the qq pair recombines to form a virtual photon. This is 
shown schematically in Fig. ^ 

The amplitude for the elastic process 7*p — > -y*p, A 1 * p {x, Q, A), is simply the product of 
amplitudes of these three subprocesses integrated over the dipole variables r and z: 

A^(x,Q,A) = J2Yl f d2r [ 1 ^n-k^z,Q)A q ,(x,r,A)^ h - h (r,z,Q), (1) 
/ h,h io 1 

where ^hh{ r i z, Q) denotes the amplitude for the incoming virtual photon to fluctuate into a 
quark-antiquark dipole with helicities h and h and flavour /. We suppress here references to 
the photon helicities for simplicity. A qq -(x, r, A) is the elementary amplitude for the scattering 
of a dipole of size r on the proton, A denotes the transverse momentum lost by the outgoing 
proton, and x is the Bjorken variable. Note that, following Ref. [1], we choose a slightly different 
convention from that commonly used, in that we include a factor of 1/ (4ir) in the integration 
measure; this convention is reflected in the normalisation of the photon and vector meson wave 
functions. 
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The elementary elastic amplitude A qg - is defined such that the elastic differential cross section 
for the qq pair scattering on the proton is 

^f = -^\A q ,(x,r,A)\ 2 , (2) 
dr 1d7t 

where t = —A 2 . It can be related to the S'-matrix element S(x,r,b) for the scattering of a 
dipole of size r at impact parameter b: 

A q g(x, r, A) = J d 2 b e- [b - A Agg{x, r,b)=i J d 2 b e~ i6A 2 [1 - S{x, r, b)] . (3) 

This corresponds to the intuitive notion of impact parameter when the dipole size is small 
compared to the size of the proton. The optical theorem then connects the total cross section for 
the qq pair scattering on the proton to the imaginary part of the forward scattering amplitude: 



(Tqq[X, r 



) = I m ^, r ,A = 0) = /d^2 [ l-Re S (,, r ,m. (4) 



The integration over b of the S-matrix element motivates the definition of the qq-p differential 
cross section as 

^M = 2[l-ReS(x,r,b)}- (5) 

The total cross section for 7*p scattering, or equivalently F 2 , is obtained, using (JTJ) and (J3J), 
by integrating the dipole cross section with the photon wave functions: 

<r0s,Q) = lmAf p L (x,Q,A = 0) = J2j d * r ^(Wfc^r), (6) 
with the overlap of the photon wave functions (\I/*\I/)£ L defined as 

= I [^k,A= + l*^A=+l + *k,A=-l*^,A=-l] , (7) 

h,h 

= j2 n hx = * h n,x=o, (8) 



h.h 



where A denotes the photon helicity and / the flavour of the qq pair. The dependence on the 
quark flavour / is specified below in Sect. 12.11 In the perturbative region, that is, for small 
dipole sizes r, the dipole cross section corresponds to exchange of a gluon ladder; see Fig. 
(left). The same diagram applies for exclusive final state production if the wave function of 
the outgoing virtual photon is replaced by the wave function of a specific final state; see Fig. El 
(right). 

The amplitude for production of an exclusive final state E, such as a vector meson (E = V) 
or a real photon in DVCS (E — 7), is given by 

Aiy Ep (x,Q,A) = Jd 2 rJ^(^)T, L A qq (x,r,A) (9) 
= i |d 2 r j^jd 2 b^* E ^ L e- ib ^2[l-S(x,r,b)], (10) 




Figure 2: The elastic scattering amplitude for inclusive DIS (left) and vector meson production 
(right). For DVCS, the outgoing vector meson in the right-hand diagram is replaced by a real 
photon. 



where (^%^)t,l denotes the overlap of the photon and exclusive final state wave functions. For 
DVCS, the amplitude involves a sum over quark flavours. This expression, used in the analysis 
of exclusive J/ip photoproduction by Kowalski and Teaney [1], is derived under the assumption 
that the size of the quark-antiquark pair is much smaller than the size of the proton. The 
explicit perturbative QCD calculation of Bartels, Golec-Biernat and Peters [40] shows that 
the non-forward wave functions can be written as the usual forward wave functions multiplied 
by exponential factors exp[±i(l — z)r ■ A/2]. Effectively, the momentum transfer A should 
conjugate to b + (1 — z)r, the transverse distance from the centre of the proton to one of the 
two quarks of the dipole, rather than to b, the transverse distance from the centre of the proton 
to the centre-of-mass of the quark dipole; see the right-hand diagram of Fig. |21 

Assuming that the S-matrix element is predominantly real we may substitute 2[1 — S(x, r, b)] 
in (Unj) with da q q/d 2 b. 

These two changes lead to 



A£r Ep (x, Q, A) = i / d 2 r f 1 ^- f d 2 6 {^)t,l e-^i-^-A ^ 



(11) 



The elastic diffractive cross section is then given by 



dt 



16n 



A. 



7*p— >Ep 
T.L 



167T 



Arc 



d 2 b 



■ (12) 



This is the basic equation for the simultaneous analysis of different exclusive processes per- 
formed in this paper. 



2.1 Forward photon wave functions 

The forward photon wave functions were perturbatively calculated in QCD by many authors; 
see, for example, Refs. [5,41]. The normalised photon wave function for the longitudinal photon 
polarisation (A = 0) is given by [9] 

* hh , x=0 (r,z,Q) = e f e^/W c 5 K _- h 2Qz{l - z)^^-, (13) 
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and for the transverse photon polarisations (A = ±1) by 



* h h,x=±i{r,z,Q) = ±e f e x/W c {ie ±Wr {z6 h: ±6- h , T - (1 - z)5 h ^8 h± ]d r + m f 5 h ,±5 h± } ~^- L , 

(14) 

where e = y/Aira cm , the subscripts h and h are the helicities of the quark and the antiquark 
respectively and 9 r is the azimuthal angle between the vector r and the x-axis in the transverse 
plane. K is a modified Bessel function of the second kind, e 2 = z(l — z)Q 2 + m 2 and N c = 3 is 
the number of colours. The flavour / dependence enters through the values of the quark charge 
e/ and mass rrif, and d r K (er) = — efTi(er). 

2.1.1 Total DIS cross sections 

In the case of the total DIS cross section cr 7 * p , which is obtained from the elastic 7*p — > 7*p 
amplitude via the optical theorem, the squared photon wave functions summed over the quark 
helicities for a given photon polarisation and quark flavour are given by the tree-level QED 
expressions: 

zf]e 2 Kl{er)+m)Kl{er)}, (15) 

z) 2 K 2 (er). (16) 

At small dipole sizes these expressions are well motivated since they can be derived from the 
LO fc r factorisation formulae. At large dipole sizes the wave functions are suppressed, since 
for large values of the argument the modified Bessel functions behave as K (er), Ki(er) ~ 
(2er) exp(— er). At larger Q 2 values the wave functions are suppressed for large r unless 
z is close to the end-point values of zero or one. 1 Near the end-points or at small Q 2 the wave 
functions are sensitive to the non-zero quark masses rrif, which prevent the integrals over r of 
the modified Bessel functions from diverging. Of course, near the end-points or at small Q 2 the 
expressions (fi3j) and (fTB^) should be considered as a model in which the value of the light quark 
masses provides a cut-off scale which should be related to the physical cut-off scale generated 
by confinement effects. It is therefore customary in dipole models to identify the light quark 
masses with the pion mass. 

2.1.2 Deeply virtual Compton scattering 

In addition to the total DIS cross section a 1 p , the photon wave functions determine also 
the DVCS process, 7*p — > 'jp. Here the outgoing photon is real and therefore the process is 
directly observed at HERA. For real photons, only the transversely polarised overlap function 

^^This is the origin of the statement that the transverse cross section is more inherently non-perturbative 
than the longitudinal cross section, since the contribution from the end-points is suppressed for the longitudinal 
but not the transverse case, see (|T5)l and (|TT)|) . 
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contributes to the cross section. Summed over the quark helicities, for a given quark flavour / 
it is given by 

2N 

= -^a em e 2 { [z 2 + (1 - zf] eK^nifK^nifr) + m 2 K (er)K (m f r)} . (17) 
2.2 Forward vector meson wave functions 

Various conventions are used in the literature for the forward vector meson wave functions. 
Recently, Forshaw, Sandapen and Shaw (FSS) [9] suggested some guidelines for bringing order 
into this problem. We will adopt their prescription in this section, apart from the overall 
normalisation factor of 1/(4-71") discussed previously, which in our case appears in the integration 
measure. 

The simplest approach to modelling the vector meson wave function is to assume, following 
Refs. [1,5,9], that the vector meson is predominantly a quark-ant iquark state and that the 
spin and polarisation structure is the same as in the photon case. In complete analogy to the 
transversely polarised photon wave function (JHJ), the transversely polarised vector meson wave 
function is 

^lh,x=±i^ z ) = ± ^^ JrT^z) { ie±Wr i z6 h,±k, T - (1 - z)d h ^5- ht± ]d r + m f 5 h> ±5- hj± } <j) T {r,z). 

(18) 

The longitudinally polarised wave function is slightly more complicated due to the fact that 
the coupling of the quarks to the meson is non-local, contrary to the photon case [9]. It is given 

by 



My + S f 



m 2 f - V 2 



Mr,z), (19) 



M v z(l - z) 

where V 2 = (l/r)d r + d 2 and My is the meson mass. The difference in the structure of 
the longitudinal wave function is due to the non-local term proportional to 5, which was first 
introduced by Nemchik, Nikolaev, Predazzi and Zakharov (NNPZ) [2,4]. 

Formulae ()18|) and ()19|) uniquely define the scalar part of the vector meson wave function 
4>T,L{ r , z)i which is obtained from the photon wave function by the replacement 

e f ez(l-z)^^ — 0r,L(r,z), (20) 

with the prefactor 2Q — > My for the case of the longitudinal polarisation. Note that this 
definition of 0t,l(^, z)|r=o matches, up to a constant factor, the definition of the distribution 
amplitude in QCD. 

The overlaps between the photon and the vector meson wave functions read then: 
(*Wt = e f e ,f c {m 2 f K Q (er)Mr,z) ~ [z 2 + (1 - z) 2 ] eK x (er)d r <fr(T,z)} , (21) 

71Z[ 1 — Z 



N r . . rn 2 - V? 



(V* V V) L = e f e — 2Qz(l - z) K (er 

7T 



(22) 
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where the effective charge if = 2/3, 1/3, or l/y/2, for J/ip, <P, or p mesons respectively. 
Although it seems to be more natural to set 5 = 1 as it was done in Refs. [2,4,9], we shall also 
use the value 8 = in order to match the assumptions of other models [1,5]. Note that the 
additional factor of 1/ [z (1 — z)\ in (|2*Tj) and (|2*2*|) as compared to the photon overlap functions 
(USD and JIHD is due to the identification (JUj). 

The usual assumption that the quantum numbers of the meson are saturated by the quark- 
antiquark pair, that is, that the possible contributions of gluon or sea-quark states to the wave 
function may be neglected, allows the normalisation of the vector meson wave functions to 
unity: 

- 1 dz 



h,h 



d 2 r 



Air 



(23) 



Thus, in the scheme presented here the normalisation conditions for the scalar parts of the 
wave functions are 



1 



Nc 
2tt 

1 



dz 



o z\\-z, 

Nc f 1 
2Wo 



J d 2 r {m}<& + [z 2 + (1 - z) 2 ] (d^r) 2 } , 
m 2 f - V 2 



dz d z r 



Mi 



V9L 



+ S 



f 



M v z(l - z) 



(24) 
(25) 



Another important constraint on the vector meson wave functions is obtained from the decay 
width. It is commonly assumed that the decay width can be described in a factorised way; 
the perturbative matrix element qq — > 7* — > factorises out from the details of the wave 
function, which contributes only through its properties at the origin. 2 The decay widths are 
then given by 



fv,T - e f 



dz 



2ttM v J z 2 (1 



V,L 



71 



dz 



{m} - [z 2 + (1 - z) 2 ] V r 2 }0 T ( 



r,z 



r=0 



M v + 5 



m. 



M v z{l - z)_ 



<f>i(r,z) 



(26) 
(27) 



r=0 



The coupling of the meson to the electromagnetic current, fv, 
electronic decay width by 

y 2 f- 



Ana 2 * 2 



3M 



is obtained from the measured 



(28) 



v 



In order to complete the model of the vector meson wave function the scalar parts of the 
wave functions 4>T,L( r , z ) should be specified. In the photon case the scalar part is given by 
modified Bessel functions, whereas for vector mesons various quark models tell us that a hadron 
at rest can be modelled by Gaussian fluctuations in transverse separation. The proton wave 
function is also directly seen to have a Gaussian form from the t-distributions of vector mesons 
at HERA; see the discussion of the proton shape below. After assuming a Gaussian form the 
modelling freedom reduces to the choice of a fluctuating variable. 



2 Usually, one assumes that the factorisation holds and that the perturbative QCD corrections are similar 
for the process of vector meson production 7*(Q 2 ) + 2g — > V and for the vector meson decay V — > 7* — * l + l~ , 
thus the corrections can be absorbed into the wave function. 
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Meson 


My/GeV 


fv 


m f /GeV 


N T 


R 2 T /GeV- 2 




R 2 L /GeV- 2 




3.097 


0.274 


1.4 


1.23 


6.5 


0.83 


3.0 




1.019 


0.076 


0.14 


4.75 


16.0 


1.41 


9.7 


p 


0.776 


0.156 


0.14 


4.47 


21.9 


1.79 


10.4 



Table 1: Parameters of the "Gaus-LC" vector meson wave functions. 



Dosch, Gousset, Kulzinger and Pirner (DGKP) [5] made the simplest assumption that the 
longitudinal momentum fraction z fluctuates independently of the transverse quark momentum 
k, where k is the Fourier conjugate variable to the dipole vector r. In what follows, this type 
of scalar wave function will be called the factorised wave function. In the DGKP model the 
parameter 5 = in (|22|h (|25|) and (|27|). The DGKP model was further simplified by Kowalski 
and Teaney [1], who assumed that the z dependence of the wave function for the longitudinally 
polarised meson is given by the short-distance limit of z(l — z) [17]. For the transversely 
polarised meson they set 4>t{t 1 z) oc [z(1 — z)] 2 in order to suppress the contribution from the 
end-points (z — > 0, 1). This leads to the "Gaus-LC" [1] wave functions given by 3 

<fa(r,z) = N T [z(l - z)] 2 exp(-r 2 /2i^), (29) 
<j> L {r, z) = N L z{l - z) ex P (-r 2 /2R 2 L ). (30) 

The values of the constants Nt,l and Rt,l m (|2*9")l and (|3*0|l . determined by requiring the correct 
normalisation and by the condition f v = f VT = f VtL , are given in Table [TJ 

The main advantage of the factorised wave functions is their simplicity. Probably a more 
realistic approach starts from the observation of Brodsky, Huang and Lepage [42] that the 
fluctuation of the quark three-momentum p in the rest frame of the meson could be described 
in a boost-invariant form. In the meson rest frame, the momentum p is connected to the qq 
invariant mass by M 2 = 4(p 2 + m 2 ). In the light-cone frame, the qq invariant mass is given by 
M 2 = (k 2 + m))/[z{\ - z)}. This leads to 

k 2 + m 2 

P= Az^z)' m ^ (31) 
and a simple ansatz for the scalar wave function in momentum space of 



4> T ,L(k,z) oc exp 



K 2 fk 2 + m 2 f 

— J - - Ami 

zl-z) f 



(32) 



This is the basis for the "boosted Gaussian" wave function of FSS [9], which was first proposed 
by NNPZ [2, 4]. 4 In the configuration space these wave functions are given by the Fourier 
transform of (f3"2~|) : 

. . - . . / m 2 n 2 2z(l-z)r 2 m 2 n 2 \ 
Mr, z) = Af TtL z(l ~ z) exp (- ^J-_ - + ^_ j . (33) 



3 Kowalski and Teaney [1] used a somewhat different convention; see the appendix for more details. 
4 Following FSS [9] we set the Coulombic part of the NNPZ wave function [2,4] to zero to avoid singular 
behaviour at the origin. This should be reasonable for p and <j> mesons, but has less justification for J ftp mesons. 
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Meson 


My/GeV 


fv 


m//GeV 


Mr 


Ml 


^ 2 /GeV- 2 


fv,T 




3.097 


0.274 


1.4 


0.578 


0.575 


2.3 


0.307 




1.019 


0.076 


0.14 


0.919 


0.825 


11.2 


0.075 


p 


0.776 


0.156 


0.14 


0.911 


0.853 


12.9 


0.182 



Table 2: Parameters of the "boosted Gaussian" vector meson wave functions. 



Note that the "boosted Gaussian" wave function has the proper short-distance limit, ~ z(l—z), 
for rrif — > 0. Following the authors of the model we set 5 = 1 in equations (|22|). ()25|) and 
(|2*7jl . defining the longitudinally polarised overlap, the normalisation and the decay constant 
respectively. We choose the "radius" parameter 7Z to reproduce the experimentally measured 
leptonic decay width of the vector meson for the longitudinally polarised case. This means that 
the calculated decay width for the transversely polarised case will be slightly different. The 
parameters 71 and Mt,l are determined by the normalisation conditions (|2*%j) and (|2*5j) and the 
decay width condition (|2~7jl . 

The parameters of the "boosted Gaussian" wave function are given in Table |2l where we 
also show the value of fv t r (EE!) computed using the given values of 7Z and Mr- (Recall that 
we require that fv,L — fv-) 

The "boosted Gaussian" wave function is very similar to the "Gaus-RF" wave function used 
in the KT investigation [1], except for the Jacobian of the transformation from the rest frame 
variables to the light-cone variables. We focus here on the "boosted Gaussian" version because 
of the proper short distance limit of the z dependence. The "CORNELL" wave function used in 
Ref. [1] cannot be used for light vector mesons since it was obtained within the nonrelativistic 
bound-state model. 

Comparing the values of the radius parameters given in Tables Q and 121 we note that the 
meson description with the "boosted Gaussian" wave function is more self-consistent; the values 
of the radius parameters Rt and Rl for the "Gaus-LC" wave functions are very different 
indicating that there are large dynamical corrections to at least one of the meson polarisation 
states. For the "boosted Gaussian" there is only one radius parameter 7Z, since the description 
of the meson is assumed to be boost-invariant between the meson rest frame and the light-cone 
frame. The shortcoming of this approach is that the predicted decay constant fv differs slightly 
between the transverse and the longitudinal polarisation components. However, the differences 
between the decay constants of the "boosted Gaussian" wave function are relatively small 
compared to the differences between the radii of the "Gaus-LC" wave function. To quantify 
this effect we fix the parameter Rt of the "Gaus-LC" wave function to the same value as Rl, 
then we predict the value of the decay constant fy,r (allowing for Nr to be determined from 
the normalisation constraint). The resulting values of fv,r were 0.44, 0.13 and 0.33 for J/ip, <p 
and p mesons respectively, to be compared with the experimental values of fy (= fv,L) of 0.27, 
0.08 and 0.16. That is, the differences between f V)T and fv,L for the "Gaus-LC" wave function 
are much larger than the equivalent differences for the "boosted Gaussian" wave function; see 
Table H 

The agreement between the decay constants for the longitudinal and transverse polarisation 
with the "boosted Gaussian" wave function is particularly good for the (f> meson wave function. 
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We note, en passant, that the difference between the two decay constants /V,r and fv,L depends 
on the assumed quark mass; for the (j) meson the difference is minimal for the strange quark 
mass of 0.14 GeV, for the J/i/j meson it is minimal for the charm quark mass of 1.15 GeV, and 
for the p meson it decreases slightly with decreasing quark mass but there is still a significant 
difference even when the quark mass is set to zero. 5 

In Fig.|3]we show the overlap functions between the photon and vector meson wave functions 
integrated over z for the three different vector mesons at Q 2 values representative of the data 
discussed later in Sect. El To be precise, we plot the quantity 



The plots show that the longitudinal overlap functions for the "Gaus-LC" and "boosted Gaus- 
sian" cases are more similar than the transverse overlap functions for all three vector mesons. 
For the (f) meson there is also a good agreement for the transverse overlap function. This indi- 
cates that observable quantities for mesons computed with either the "Gaus-LC" or "boosted 
Gaussian" wave functions should be very similar, in spite of the sizable disagreement between 
i2y and i?| for the "Gaus-LC" wave function. 

2.3 Dipole cross sections 

2.3.1 Review of dipole cross sections 

The dipole model became an important tool in investigations of deep-inelastic scattering due 
to the initial observation of Golec-Biernat and Wiisthoff (GBW) [35, 36] that a simple ansatz 
for the dipole cross section integrated over the impact parameter b, a q q, was able to describe 
simultaneously the total inclusive and diffractive DIS cross sections: 



where <t is a constant and Q s (x) denotes the x dependent saturation scale, Q 2 (x) = (x /x) Xgbw 
GeV 2 . The parameters <Jq = 23 mb, Agbw — 0.29 and x$ = 3 x 10~ 4 were determined from a fit 
to the F2 data without including charm quarks. After inclusion of the charm quark contribution 
with mass m c = 1.5 GeV into the fit, the parameters of the GBW model changed to o"o = 29 
mb, Agbw = 0.28 and xo = 4 x 10~ 5 . Although the dipole model is theoretically well justified 
for small-size dipoles only, the GBW model provided a good description of data from medium 
Q 2 values (~ 30 GeV 2 ) down to low Q 2 (~ 0.1 GeV 2 ). The saturation scale Q 2 is intimately 
related to the gluon density in the transverse plane. The exponent Agbw determines therefore 
the growth of the total and diffractive cross sections with decreasing x. For dipole sizes which 
are large in comparison to 1/Q S the dipole cross section saturates by approaching a constant 
value cr , which becomes independent of Agbw- It is a characteristic feature of the GBW model 

5 For the <fi meson, the relative difference of decay constants fv,T and fv,L is 11% for m s — 0.3 GeV and 3% 
for m s = 0.05 GeV. For the p meson, the relative difference of decay constants is 36% for m u ,d — 0.3 GeV and 
14% for m u , d = 0.05 GeV. 




(34) 




(35) 
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Figure 3: Overlap functions (j2TJ and (j22| between the photon and vector meson wave functions 
integrated over z for the three different vector mesons at Q 2 values representative of the data. 
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that a good description of data is due to large saturation effects, that is, the strong growth due 
to the factor x _Agbw is, for large dipoles, significantly flattened by the exponentiation in (|35|). 

The assumption of dipole saturation provided an attractive theoretical background for in- 
vestigation of the transition from the perturbative to non-perturbative regimes in the HERA 
data. Despite the appealing simplicity and success of the GBW model it suffers from clear 
shortcomings. In particular it does not include scaling violations, that is, at large Q 2 it does 
not match with QCD (DGLAP) evolution. Therefore, Bartels, Golec-Biernat and Kowalski 
(BGBK) [32] proposed a modification of the original ansatz of (|33j) by replacing Q 2 by a gluon 
density with explicit DGLAP evolution: 



The scale of the gluon density, /i 2 , was assumed to be /i 2 = C/r 2 + /i 2 ,, and the gluon density 
was evolved according to the leading-order (LO) DGLAP equation without quarks. 

The BGBK form of the dipole cross section led to significantly better fits to the HERA F 2 
data than the original GBW model, especially in the region of larger Q 2 . The good agreement 
of the original model with the DIS diffractive HERA data was also preserved. However, the 
contribution from charm quarks was omitted in the BGBK analysis. 

The BGBK analysis found, surprisingly, that there exist two distinct solutions, both giving 
a very good description of the HERA data, depending on the quark mass in the photon wave 
function. The first solution was obtained assuming m u ^ s = 0.14 GeV and led to the initial 
gluon density, xg(x, /Zq) oc x _Afl , with the value of exponent X g = 0.28 at /ig = 0.52 GeV 2 , which 
is very similar to the Agbw- As in the original model, the good agreement with data was due 
to substantial saturation effects. In the second solution, which took m u ^ s = 0, the value of the 
exponent was very different, X g = —0.41 at a fixed /ig = 1 GeV 2 . The initial gluon density no 
longer rose at small x; it was valence-like, and QCD evolution played a much more significant 
role than in the solution with m Ut d tS = 0.14 GeV. 

The DGLAP evolution, which is generally used in the analysis of HERA data, may not be 
appropriate when x approaches the saturation region. Therefore, Iancu, Itakura and Munier [37] 
proposed a new saturation model, the Colour Glass Condensate (CGC) model, in which gluon 
saturation effects are incorporated via an approximate solution of the Balitsky-Kovchegov 
equation [43-45]. The CGC dipole cross section is 



where Y = ln(l/x), 7 S = 0.63, k = 9.9 and Q s = Q s (x) = (x /x) x ^ 2 . The free parameters a , 
Ao, A and Xq were determined by a fit to HERA F 2 data. The coefficients A and B in the 
second line of (|3T|) are determined uniquely from the condition that cr^ 30 , and its derivative 
with respect to rQ s , are continuous at rQ s = 2: 



BGBK 



(x,r) = (T {l - exp [-7i 2 r 2 a s ([i 2 )xg(x, fi 2 ) / (3cr )] } . 



(36) 




(37) 



.4 




1 , (i-A/p 



(38) 



(l-jV ) 2 ln(l-jV ) 
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Later, also Forshaw and Shaw (FS) [38] proposed a Regge-type model with saturation effects. 
The CGC and FS models provide a description of HERA inclusive and diffractive DIS data 
which is better than the original GBW model and comparable in quality to the BGBK analysis. 
Both models find strong saturation effects in HERA data comparable to the GBW model and 
the solution of the BGBK model with m Uj d, s = 0.14 GeV. 

All approaches to dipole saturation discussed so far ignored a possible impact parameter 
dependence of the dipole cross section. This dependence was introduced in this context by 
KT [1], who assumed that the dipole cross section is a function of the opacity Q, following for 
instance Ref. [3]: 

^ = 2 (!-.-»). (39) 

At small x the opacity Q can be directly related to the gluon density, xg(x,fi 2 ), and the 
transverse profile of the proton, T(b): 



TT 



2 



n = —r 2 a s (^)xg(x^ 2 )T(b). (40) 



N, 



c 



The formulae of (|39|) and (J4(Jj) are called the Glauber-Mueller dipole cross section. The diffrac- 
tive cross section of this type was used around 50 years ago to study the diffractive dissociation 
of deuterons by Glauber [46] and reintroduced by Mueller [47] to describe dipole scattering in 
deep-inelastic processes. 



2.3.2 Applied dipole cross sections 



Since the description of exclusive vector meson production is the focus of this investigation we 
concentrate here on impact parameter dependent dipole cross sections. First, we use the same 
form of the differential dipole cross section as in the KT investigation [1]: 



da. 



gq 



d 2 b 



1 — cxp 



7T 

"2N, 



■r 2 a s (n 2 )xg(x,n 2 )T(b) 



(41) 



Here, the scale /u, 2 is related to the dipole size r by /i 2 = 4/r 2 + /x 2 ,. The gluon density, xg(x, (i 2 ), 
is evolved from a scale /i 2 , up to /i 2 using LO DGLAP evolution without quarks: 

dxg{x,jj?) ois{n 2 ) 



/ dzP gg {z)-g{- ^ 

J x Z \Z 



<91n/i 2 2n 
The initial gluon density at the scale /i 2 , is taken in the form 

xg(x,fj%) = A g x- Xs (1 -x) 5 - 6 



(42) 



(43) 



The values of the parameters /J 2 ,, A g , and X g are determined from a fit to Fi data. For the 
light quarks, the gluon density is evaluated at x = x B (Bjorken-x), while for charm quarks, 
x = x B (l + 4m 2 /Q 2 ). The contribution from beauty quarks is neglected. For vector meson 
production, the gluon density is evaluated at x = x B (l + My/Q 2 ). The LO formula for the 
running strong coupling as(/i 2 ) is used, with three fixed flavours and Aqc D = 0.2 GeV. 
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The proton shape function T(b) is normalised so that 

d 2 b T(b) = I. (44) 



We consider first a Gaussian form for T(b), that is, 

T G (b) = ai-e-^fe, (45) 

where Bq is a free parameter which is fixed by the fit to the differential cross sections da/dt 
for exclusive vector meson production. This distribution yields the average squared transverse 
radius of the proton, 

(b 2 ) = 2B G . (46) 

Assuming that the Gaussian distribution given by (jUj) holds also in three dimensions (with a 
different normalisation factor) we obtain the relationship between the parameter Bq and the 
Hofstadter radius of the proton R p , namely R 2 = 3Bq- Note that the Hofstadter experiment [48] 
measured the electromagnetic radius whereas we probe the gluonic distribution of the proton. 

The two-dimensional Fourier transform of (J43J) has the exponential form which is supported 
by the data: 6 

oce-^W. (47) 

dt v ; 

Alternatively, we assume that the gluonic density in the proton is evenly distributed over a 
certain area within a sharp boundary, and is zero beyond this boundary. That is, we assume a 
step function, again normalised as in 



T siP) = (b s - b) , (48) 

where b$ is a free parameter, for which the average squared transverse radius of the proton is 

(b 2 ) = |. (49) 

This is the form of T(b) implicitly used in all b- independent parameterisations of the dipole 
cross section. That is, it is usually assumed that 

^ = 2[1 -BeS(x,r,b)] = 2 Af(x, r,b) = 2 Af(x, r) 6 (b s - b) , (50) 
so that integration over b gives 



(x,r) = a Af(x,r), (51) 



where the parameter o"o = 2irb 2 s is usually obtained by fitting to the F2 data. This is the form 
assumed in the GBW model (J^J), the BGBK model flU, and the CGC model fl37|). Note that 

6 Note that for exclusive diffractive processes at large values of (My + Q 2 ) the typical dipole size r is small, 
and the ^-dependence of the cross section is determined entirely by the proton transverse profile. 
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the scattering amplitudes Af(x, r, b) or Af(x, r) can vary between zero and one, where Af = 1 is 
the unitarity limit. 

To introduce the impact parameter dependence into the CGC model [37], we modify (p?Tj) 
to obtain the "b-CGC" model: 



da, 



d 2 b 



= 2Af(x,r,b) = 2 x 



1 _ e -A\n 2 {BrQ s ) 



rQ s < 2 
rQ s > 2 



(52) 



where now the parameter Q s depends on the impact parameter: 



Q s = Q s (x,b) = ( — 
x 



exp 



2B, 



CGC 



273 



(53) 



Note that, in contrast to the parameter Bq in the KT approach, a straightforward interpretation 
of -Bccc hi terms of the proton size is not possible due to the r and Y dependence of the exponent 

2 ( 7s + ^ ln ir) in 

Following KT [1] we define the saturation scale Q 2 S = 2/r|, where the saturation radius r$ 
is the dipole size where the scattering amplitude Af has a value of 1 — exp(— 1/2) ~ 0.4, that 
is, rs is defined by solving 

Af{x,r s ,b) = l-e~3, (54) 

with the same condition for the b- independent dipole models. For the GBW model (J35|) . the 
saturation scale Q 2 S — 2/r| defined by (|H1|) coincides with Q 2 s (x) = (xo/x) Xgbw GeV 2 . However, 
for the CGC (13*7)) and b-CGC models, the saturation scale Qs defined by (jo^j) differs from 
the parameter Q s . Note that we use upper-case S and lower-case s to distinguish between these 
two scales. The saturation scale Qs is the quantity we shall later compute and compare for the 
different dipole models in Sect. El 



2.3.3 Phenomenological corrections for exclusive processes 

After performing the angular integrations, (jllj) reduces to 

/OO pi J /"OO J 

dr (2?rr) ^ ^ ^ d6 (27Tb) (n*h,L MbA) J ([l - «]rA) (55) 

where Jo is the Bessel function of the first kind and E — V, 7 denotes either the exclusive vector 
meson or DVCS final state. The derivation of the expression for the exclusive vector meson 
production or DVCS amplitude, (JTTJ), relies on the assumption that the S'-matrix is purely real 
and therefore the exclusive amplitude A is purely imaginary. The real part of the amplitude 
can be accounted for by multiplying the differential cross section for vector meson production 
or DVCS, (H2J), by a factor (1 + (3 2 ), where (3 is the ratio of real to imaginary parts of the 
scattering amplitude A and is calculated using 

d In ( X / * P ^ E A 

(3 = tan(7rA/2), with A = \ —. J . (56) 

a mil x) 
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This procedure (or similar) is adopted in other descriptions of vector meson production to 
account for the real part of the amplitude; see, for example, Refs. [4,9,19]. 

For vector meson production or DVCS, we should use the off-diagonal (or generalised) gluon 
distribution, since the two gluons in the right-hand diagram of Fig. |21 carry different fractions 
x and x' of the proton's (light-cone) momentum. In the leading ln(l/x) limit, the skewed effect 
vanishes. However, the skewed effect can be accounted for, in the limit that i' < i < 1, by 
multiplying the gluon distribution xg(x,/i 2 ) in (jUJ) by a factor R g , given by [49] 

2^+ 3 r(A + 5/2) d\n[xg(x^ 2 )] 
^ (A) ~ A T(A + 4) ' Wlth A= dln(l/x) ■ (57) 

This skewedness effect is also accounted for in the calculation of vector meson production 
by Martin, Ryskin and Teubner (MRT) [19], but is neglected in most other dipole model 
descriptions. 



3 Description of HERA data with the "b-Sat" model 

In this section we describe HERA data within the generalised impact parameter dipole satu- 
ration ( "b-Sat" ) model in which the dipole cross section is given by (|4Tj) and the proton shape 
function T{b) is assumed to be purely Gaussian (J43j) . The total DIS cross section is given by © 
and the photon overlap functions by (JX3J) and (fTtjj). For exclusive processes, the differential cross 
sections are given by ()12[) with the phenomenological improvements described in Sect. 12.3.31 
For vector mesons the overlaps of wave functions are given by (f2*TJ) and ()22j) . and for the DVCS 
process by (JT7J). 

The light quark masses are taken to be m Uj d, s = 0.14 GeV, the value of the pion mass, 
which ensures the proper exponential cut-off of the photon wave functions (fT5|) and (|T6|) at 
large distances. The value of the charm mass was chosen to be m c = 1.4 GeV, but other 
choices for the charm and light quark masses are also discussed below. The free parameters of 
the model are /1q, A g and X g of the initial gluon distribution, xg(x, /1q) = A g x~ Xg (1 — x) 5,6 , and 
the proton width Bq. The aim of the model is to describe with these four parameters the total 
DIS cross section for x B < 0.01 and all total and differential cross sections for J/ip, <fi and p 
meson production, as well as DVCS. The dipole cross section as determined in the b-Sat model 
is shown at various impact parameters in Fig. 



3.1 Total 7*p cross section 

The parameters in the initial gluon distribution fi3|) are determined by fitting the ZEUS F 2 
data [22,23] with x B < 0.01 and Q 2 e [0.25,650] GeV 2 . They are obtained in a quickly 
converging iterative procedure in which the F 2 data are fitted alternately with the t-distributions 
of the vector meson data (see Sect. 13. 3J) which determine the parameter Bq = 4 GeV -2 . As 
well as our main fit with m u ^ )S = 0.14 GeV and m c = 1.4 GeV, shown in the first line of Table 
I2J we also make alternative fits with different quark masses. As in Ref . [1] , the best fit to F 2 is 
obtained with very low light quark masses, m u ^ iS — 0.05 GeV. The quark mass of 0.14 GeV, 
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Figure 4: Dipole cross section at various impact parameters, as determined in the b-Sat model. 



Model 


T(b) 


Q 2 /GeV 2 


m uAs /GeV 


m c /GeV 


H 2 /GeV 2 


A 9 




X 2 /d.o.f. 


b-Sat 


Gaussian 


[0.25,650] 


0.14 


1.4 


1.17 


2.55 


0.020 


193.0/160 = 1.21 


b-Sat 


Gaussian 


[0.25,650] 


0.14 


1.35 


1.20 


2.51 


0.024 


190.2/160 = 1.19 


b-Sat 


Gaussian 


[0.25,650] 


0.14 


1.5 


1.11 


2.64 


0.011 


198.1/160 = 1.24 


b-Sat 


Gaussian 


[0.25,650] 


0.05 


1.4 


0.77 


3.61 


-0.118 


144.7/160 = 0.90 


b-Sat 


Step 


[0.25,650] 


0.14 


1.4 


1.50 


2.20 


0.071 


199.6/160 = 1.25 



Table 3: Parameters of the initial gluon distribution ()43|) determined from fits to F 2 data [22,23]. 
All predictions using the b-Sat model in this paper are evaluated with the set of parameters 
given in the first line unless explicitly stated otherwise. 
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Figure 5: Top: The total DIS cross section o"7 ot P vs. iy 2 for different Q 2 . The data points 
plotted are from ZEUS [22,23]. Bottom: The A to t parameter for inclusive DIS defined by 
a£t oc (l/x) Atot . The data points plotted are from ZEUS [22,23] and HI [24]. 
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Figure 6: The gluon distribution xg(x,fi 2 ) for different fixed p 2 (left) or fixed r (right). 



which is more appropriate as a cut-off mass for vector meson bound states, gives a fit to F 2 of 
still acceptable quality; see Table El The last line of Table El shows also the fit results performed 
with the step-like proton shape defined by (J48j) with the parameter 65 = 4 GeV -1 , which we 
discuss further in Sect. El 

To compare with the fits obtained by global analysis using the NLO DGLAP formalism, 
we evaluated the \ 2 for a subset of the ZEUS F 2 data [23] with x B < O.Of and Q 2 > 2 GeV 2 
comprising of 116 data points. The main b-Sat fit shown in the first line of Table El gave a x 2 
of 114, while the most recent NLO DGLAP fit by the MRST group [50] gave a \ 2 of 96 for the 
116 data points. 

In Fig. El we show the comparison of the main b-Sat fit results with measurements of the 
total DIS cross section <7^ ot p . In the same figure we also show the comparison for the rate of 
rise of the total DIS cross section, A tot , defined by oc (I/x) Atot . Both comparisons show a 
very good agreement between data and the b-Sat model results. 

Let us make some general remarks about the sensitivity of the fit to the assumed quark 
masses and the proton shape. Table El shows that the variation of the charm quark mass 
does not sizably change the fit parameters. On the other hand, the choice of the light quark 
mass influences the value of the X g parameter and consequently the evolution of the gluon 
density. In Fig. El we show the gluon distribution for different scales p 2 or dipole sizes r. The 
correlation between the assumed value of the light quark mass and the X g and /j,q parameters 
was investigated in detail in Ref. [1]. Consequently, in the b-Sat model the description of the 
change of the parameter A to t with Q 2 is mainly due to evolution effects and not to saturation 
effects as in, for example, the GBW model [35,36]. 

3.2 Vector meson total cross sections 

We now compare our predictions for exclusive vector meson production with recent published 
HERA data for J ftp [25-27], <ft [28] and p [29] meson production. 7 The HI J/ip cross sec- 

7 The ZEUS j*p cross sections [25,26,28] are given as a = ot + &l, while HI [27,29] give a = <jt + £&l, 
where e = (1 — y)/(l — y + J/ 2 /2) and (e) as 0.99. We use the ZEUS definition in our calculations. 
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Figure 7: Total vector meson cross section a vs. (Q 2 + M v ) compared to predictions from the 
b-Sat model using two different vector meson wave functions. The ZEUS J / ip photoproduction 
point is taken from Table 1 of Ref. [25], from the muon decay channel with W = 90-110 GeV. 

tions [27] are measured in the range \t\ < 1.2 GeV 2 while ZEUS measure \t\ < 1 GeV 2 for 
electroproduction [26] and \t\ < 1.8 GeV 2 (J/ip -> p+ir) or \t\ < 1.25 GeV 2 (J/ip -> e+e - ) 
for photoproduction [25]. The ZEUS <p data [28] have \t\ < 0.6 GeV 2 , while the HI p data [29] 
have |t| < 0.5 GeV 2 . 

In Fig. [7| we show the (Q 2 + My) dependence of the total cross section a for all three 
vector mesons at a fixed value of W. The inner error bars indicate the statistical uncertainties 
only, while the outer error bars include the systematic uncertainties added in quadrature. The 
predictions are given integrated over the appropriate t range. For the J/ip data, the predictions 
shown correspond to the HI t range. The predictions of the model are in good agreement 
with data for both vector meson wave functions. The model reproduces the Q 2 dependence as 
well as the absolute magnitude of the data. The prediction for the absolute normalisation is 
determined mainly by the gluon density obtained from the fit to the total DIS cross section 
(or F 2 ) and the shapes of the "Gaus-LC" and "boosted Gaussian" wave functions, discussed 
in Sect. I2.21 Although these two vector meson wave functions are quite different, they lead 
to similar predictions using the constraints from the normalisation and vector meson decay 
width conditions given in (|2^ . (|55|) . and (|27jl . Note that, unlike the MRT calculations [19] 
compared to the HI J/ip data in [27], we do not require an additional normalisation factor 
~ 2 to achieve agreement with the data. Note also that the MRT calculations [19], based 
on ^-factorisation using an unintegrated gluon distribution, take as input the gluon density 
determined from the global analyses using collinear factorisation. There is no a priori reason 
why the fitted parameters in the two gluon distributions determined in these two calculational 
frameworks should be identical. The dipole approach is self-consistent in that the gluon density 
is determined from the inclusive process and applied to exclusive processes within the same 
calculational framework. 
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Figure 8: Total vector meson cross section a vs. W compared to predictions from the b-Sat 
model using two different vector meson wave functions. The ZEUS J/ip data points [25,26] 
have been scaled to the HI Q 2 values [27] using the Q 2 dependence measured by ZEUS of the 
formaoc (Q 2 + M 2 )" 2 - 44 [26]. 

In Fig. Owe show the W dependence of the total cross section a for fixed values of Q 2 . 
Here, the "boosted Gaussian" vector meson wave function gives a slightly better description 
of the data. In Fig. |U] we show the effect of changing the charm quark mass from the default 
value of 1.4 GeV to 1.35 GeV or 1.5 GeV. We also show the effect of changing the light quark 
masses from 0.14 GeV to 0.05 GeV. In each case, we refit the Fi data to determine the gluon 
distribution with parameters given in Tabled The absolute magnitude of the J ftp cross sections 
is strongly dependent on the choice of the charm quark mass, particularly at small Q 2 values. 
The cross sections for the and p vector mesons are only weakly dependent on the choice of 
the light quark masses. This is because, in the Q 2 range considered in this paper, the scale for 
light vector meson production, given by e 2 = z(l — z)Q 2 + m 2 , is predominantly given by Q 2 
whereas for J / ip mesons the scale e 2 is dominated by the square of the charm quark mass. Note 
also that for all vector mesons the sensitivity of the cross section to the quark mass decreases 
with increasing Q 2 . 

We then perform a fit to the theory predictions shown in Fig. |H] of the form o oc W s and 
compare the values of 5 obtained to the experimental values; see Fig. [TUl For p production, 
we instead show atp(0) calculated from 6 = 4[a F ((i)) — 1], where ap((£)) = aip(0) + a p (t), 
(t) = —1/B D , B D is the theoretical prediction (see Fig. 114]) . and a' P = 0.25 GeV -2 . We observe 
again a reasonable agreement of the model results with data. 

A variable which is more sensitive to the details of the wave function is the ratio of the 
longitudinal to the transverse cross sections, R = a L /a Tl shown in Fig. [TTJ This is due to 
the fact that the ratio <Jl/&t probes the behaviour of the transversely polarised vector meson 
wave function close to the end-points (z — ► 0, 1). At large values of Q 2 , the contributions 
from the intermediate values of z ~ 1/2 follow the simple, perturbative scaling that leads to 
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Y*p -> J/y p y*p -> <t> p Y*P -j P P 




100 200 300 50 100 150 50 100 150 
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Figure 9: Total vector meson cross section a vs. W compared to predictions from the b-Sat 
model using the "boosted Gaussian" vector meson wave function for different quark masses. 




Figure 10: The power 5 vs. (Q 2 + My), where 5 is defined by fitting a oc W 5 , compared 
to predictions from the b-Sat model using two different vector meson wave functions. For p 
mesons, we instead show cnp(O) defined in the main text; the error bars represent the statistical 
and non-correlated systematic uncertainties only [29]. 
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Figure 11: The ratio R = <Jl/o~t vs. Q 2 compared to predictions from the b-Sat model using 
two different vector meson wave functions. 

g lI g t ~ Q 2 - This simple scaling is affected by the Q 2 evolution of the anomalous dimension of 
the gluon distribution [15,19], and by the contributions from the end-points to the transverse 
cross section, which are different for the "Gaus-LC" and "boosted Gaussian" vector meson 
wave functions. Fig. El shows that the "boosted Gaussian" wave function is favoured by the 
p meson data, where the "Gaus-LC" wave function leads to a value of <Jl/&t which rises too 
rapidly with increasing Q 2 . For J/ip and mesons, both vector meson wave functions lead to a 
similar behaviour. In Fig. we show the effect of changing the quark masses when using the 
"boosted Gaussian" wave function. For p mesons, the ratio ctl/o't shows a strong dependence 
on the quark mass. A more precise analysis, which goes beyond the scope of this paper, shows 
that the ratio (JlI^t is very sensitive to the behaviour of the wave functions at the end-points 
(z-0,1). 



3.3 Vector meson ^-distributions 

The observed ^-distributions of the vector meson processes are an important source of infor- 
mation on the shape of the proton in the low-x region. Fig. EH shows the HERA data on 
t-distributions for J/ijj [25-27] and <p [28] meson production. Fig. El shows the effective slope 
of the ^-distribution, the parameter Bp, for J/ip, <fi and p [29] vector mesons as a function of 
(Q 2 + My). The parameter Bd describes the area size of the interaction region and is ob- 
tained by making a fit to the observed (or computed in the model) t-distributions of the form 
da/dt oc exp(— B D \t\). The theory predictions for B D are all obtained by making fits to da/dt 
in the range \t\ < 0.5 GeV 2 . Figs. ITH1 and ITU show that the t dependence and the (Q 2 + My) 
dependence of Bd are well described by the dipole model predictions for all three vector mesons 
whether using either the "Gaus-LC" or the "boosted Gaussian" vector meson wave functions. 
We note that this good description is obtained with only one value of the width of the proton 
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Figure 12: The ratio R = <Jl/o~t vs. Q 2 compared to predictions from the b-Sat model using 
the "boosted Gaussian" vector meson wave function for different quark masses. 
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Figure 13: Differential vector meson cross section da/dt vs. \t\ compared to predictions from 
the b-Sat model using two different vector meson wave functions. The ZEUS photoproduction 
(J/ip — > fi + data points [25] shown in the second plot show only the statistical errors 
and are for W = 90-110 GeV with the predictions calculated at W = 100 GeV. The ZEUS 
electroproduction data points [26] shown in the same plot are for W = 90 GeV. 
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Figure 14: The t-slope parameter B D vs. (Q 2 + M v ), where B D is defined by fitting da/dt oc 
exp(— B D \t\), compared to predictions from the b-Sat model using two different vector meson 
wave functions. 

shape, Bq. 

The proton shape, in the b-Sat model, is assumed to be purely Gaussian The width 

of the Gaussian, Bq, determined by optimising the agreement between the model predictions 
and data for the t-distributions of the vector mesons and their effective slopes Bp, is found 
to be Bq = 4 GeV -2 . This value is mainly determined by the ^-distributions of J/ip mesons 
measured by ZEUS [25,26] and HI [27]. We note, however, that although the values of the 
Bd parameters measured by the two experiments are in agreement within errors, the spread of 
their values is somewhat large; see the first plot of Fig. El We estimate the error on the value 
of the parameter Bq as being around 0.5 GeV -2 . 

The value of Bq = 4 GeV -2 found in this investigation is slightly smaller than in the KT [1] 
investigation where Bq = 4.25 GeV -2 was determined using only the ZEUS J/ip photoproduc- 
tion data [25]. Fig. ITU shows that the subsequent ZEUS measurements of J/ifj electroproduc- 
tion [26] exhibit higher values of Bp and therefore require a higher value of Bq. Note that 
the effect of taking the size of the vector meson into account, that is, including the BGBP [40] 
factor in (J 12)) arising from the non-forward wave functions, exp [i(l — z)r ■ A], lowers the cross 
section for non-zero t and therefore lowers the required value of Bq; recall that this factor was 
neglected by KT [1]. 

Note also that the obtained values of Bp at the same (Q 2 + M v ) are larger for light vector 
mesons than for J/ip, in accordance with the data. This occurs because the scales Q 2 and m 2 
enter the photon wave function in slightly different ways. We shall illustrate this by comparing 
J/ijj photoproduction with light vector meson electroproduction at the same value of (Q 2 + M v ), 
implying Q 2 ~ 4m 2 . The characteristic size of the scattering dipole is set by 1/e with e 2 = 
z{l — z)Q 2 + m 2 . For the photoproduction of J ftp, e has no z dependence, e 2 = m 2 . In 
contrast, for light vector mesons e 2 varies with z from Q 2 /A at z = 1/2 down to m 2 d s at z — > 
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Figure 15: The t-slope parameter Bd vs. (Q 2 + My) compared to predictions from the b- 
Sat model using the "boosted Gaussian" vector meson wave function. We show the effect of 
switching off the eikonaiisation in the dipole cross section (|41|) . and omitting the BGBP [40] 
factor, exp [i(l — z)r -A], in 

and z — > 1, so that the effective value of e 2 is significantly lower than Q 2 /4 + m 2 ds ~ m 2 . 
Therefore, for light vector meson production at Q 2 ~ 4m 2 , the typical dipole size is larger 
than for photoproduction of J/ip. This is particularly pronounced at the end-points z —>■ 
and z — > 1 for the transversely polarised light vector mesons. At sufficiently large values of 
Q 2 , however, the longitudinally polarised mesons dominate and the typical dipole size becomes 
small enough to have a negligible contribution to Bd for both light and heavy mesons. Hence, 
at large (Q 2 + My), Bd tends to a universal value determined by the proton shape alone. 

It is important to realise that the dependence of Bd on (Q 2 + M v ) observed for light vector 
mesons originates from the enlargement of the interaction area due to the dipole transverse 
extension. Recall that this effect is taken into account by the BGBP [40] prescription of the 
QCD dipole scattering at t ^ 0. It also partly arises from the saturation effects which play a 
stronger role for the larger typical dipole sizes at small (Q 2 +M V ). We investigate the interplay 
between these two mechanisms on the value of Bd in Fig.E3 We show the effect of switching off 
the eikonaiisation, that is, replacing the dipole cross section by the opacity Q (|4*0|) . We also 
show the effect of omitting the BGBP [40] factor, exp [i(l — z)r - A], in (fT2j) . Without these 
two effects, which diminish with increasing (Q 2 + My), the t-slope Bd tends to the universal 
value of Bd = Bq = 4 GeV -2 . Without the BGBP factor, the eikonaiisation has a significant 
effect for <fi and p mesons, but it is not enough to describe the Bd data points. With the BGBP 
factor, the eikonaiisation has only a small effect and the rise of Bd with decreasing (Q 2 + My) 
nicely reproduces the rise observed in the data. 

We also investigated, for completeness, the W dependence of the t-distributions. In Fig. HB] 
we show the W dependence of da/dt for fixed values of |t| and Q 2 . For each value of t, we make 
a fit of the form da/dt oc W^ av ^~^ and then plot otp(t) against |t|; see Fig. [T7| We also fit the 
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Figure 16: Differential vector meson cross section da/dt vs. W compared to predictions from 
the b-Sat model using two different vector meson wave functions. 
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Figure 17: The Pomeron trajectory otp(t) vs. where Qfp(t) is determined by fitting da/dt oc 
yy4[ap(t)-i] ^ compared to predictions from the b-Sat model using two different vector meson 
wave functions. 
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Figure 18: The t-slope parameter B D vs. W, where B D is defined by fitting da/dt oc 
exp(— Bo\t\), compared to predictions from the b-Sat model using two different vector me- 
son wave functions. 
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Figure 19: Total DVCS cross sections a vs. Q 2 (left) and a vs. W (right) compared to predictions 
from the b-Sat model. 



same data to the form da/dt oc exp(—B D \t\) for each value of W, then we plot B D against W; 
see Fig. d 



3.4 Deeply virtual Compton scattering 

We now compare to the recently published DVCS data from HI [30] and ZEUS [31]. We use 
the b-Sat model with a Gaussian T{b) and Bq = 4 GeV -2 , and quark masses m Uy d, s = 0.14 
GeV and m c = 1.4 GeV. In Fig. (left) we show the Q 2 dependence of the cross section 
integrated over \t\ up to 1 GeV 2 for W = 82 GeV compared to the HI data [30]. We also show 
the ZEUS data [31] at W = 89 GeV rescaled to W = 82 GeV using a oc W s , with S = 0.75 [31]. 
In Fig. EI] (right) we show the W dependence of the cross section integrated over \t\ up to 1 
GeV 2 for Q 2 = 8 GeV 2 compared to the HI data [30]. We also show the ZEUS data [31] at 
Q 2 = 9.6 GeV 2 rescaled to Q 2 = 8 GeV 2 using a oc Q~ 2n , with n = 1.54 [31]. Fitting the theory 
predictions to the form a oc W s gives 5 = 0.80 to be compared with the experimental value of 
0.77 ±0.23 ±0.19 [30]. We see from Fig.DJthat the Q 2 and W dependence of the DVCS data, 
as well as the absolute normalisation, are well described by the b-Sat model. 

The t-distribution is shown in Fig. 1201 for Q 2 = 8 GeV 2 and W = 82 GeV compared to the 
HI data [30]. At small t the data are well- described, while at larger t the prediction slightly 
overestimates the data, due to a t-slope which is too small. Fitting the theory prediction to 
the form da/dt oc exp(— Br>\t\) for \t\ < 0.5 GeV 2 gives Bo = 5.29 GeV -2 , to be compared 
with the experimental value of 6.02 ± 0.35 ± 0.39 GeV -2 [30]. When comparing these values 
one should bear in mind that the value of the parameter Bq = 4 GeV -2 determined from the 
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Model Q 2 /GeV 2 


m uAs /GeV m c /GeV 


Ao 


xo/10- 4 


A 


xVd.o.f. 


b-CGC [0.25,45] 


0.14 1.4 


0.417 


5.95 


0.159 


211.2/130 = 1.62 



Table 4: Parameters of the b-CGC model, (J52j) and (J53)) . determined from a fit to F 2 data 
[22,23]. 



t-distributions of the vector meson data has a possible uncertainty which could be as large as 
0.5 GeV- 2 . 

Summarising, we can see that the agreement of the predictions from the b-Sat model with 
DVCS data is remarkably good, especially if we note that the DVCS data were not used in 
fixing any parameters of the model. 

4 Impact parameter dependent CGC model 

We have seen that almost all features of the exclusive diffractive HERA processes are well 
described by the impact parameter dependent saturation ("b-Sat") model with a Gaussian 
T(b) of width Bq = 4 GeV -2 . The b-Sat model assumes the validity of DGLAP evolution 
which may not be appropriate when x approaches the saturation region. Therefore, we also 
investigated the impact parameter dependent CGC ("b-CGC") model, in which the dipole cross 
section is given by (J52J) and (J53J). In the b-CGC model the evolution effects are included via 
an approximate solution to the Balitsky-Kovchegov equation [43-45]. 

Similar to the b-Sat model, the parameter -Bcgc = 5.5 GeV~ 2 in (|53|) is determined by 
requiring a good description of the t-slopes of vector meson data, while the three parameters 
Ao, A and xo in and (|53~|l are determined by fitting the F 2 data [22,23] with x B < 0.01 
and Q 2 G [0.25,45] GeV 2 . The results of the fit are shown in Table EJ The fit to the F 2 data 
with the b-CGC model gives a sizably worse description than the b-Sat model as seen from the 
value of the x 2 /d.o.f. in Table |U and the comparison with data of the parameter A to t shown 
in the bottom plot of Fig. The significant deterioration of the fit quality is due to the fact 
that in the impact parameter dependent description, saturation effects can only be sizable in 
the core of the proton, see the discussion in Sect. El The relatively poor quality of the fit 
is the main reason why we prefer to use a DGLAP-evolved gluon density together with the 
Glauber-Mueller dipole cross section, that is, the b-Sat model. 

Although almost all features of the vector meson and DVCS data are well described by the 
b-Sat model, there is one exception, namely a' P . It is predicted to be close to zero, due to 
the assumed factorisation of T(b) from the gluon distribution xg(x,/j, 2 ), in some disagreement 
with the data; see Figs. [T7| and [THJ In the b-CGC model the W (or x) dependence is not 
factorised from the b dependence. Therefore, an appreciable a' ¥ is achievable, as shown in 
Fig.|^ Here, we use the "boosted Gaussian" vector meson wave function in both cases. In fact, 
for photoproduction, a fit to the model predictions of the form B D = B + 4a' F \n[W/ (90 GeV)] 
gives a' F = 0.075 for the b-CGC model compared to a' F = 0.004 for the b-Sat model. However, 
the value of a' ¥ from the b-CGC model is still slightly low when compared to the values of 
0.116 ± 0.026±8$g [25] or a' P = 0.164 ± 0.028 ± 0.030 [27] measured by experiment. We note 
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Figure 21: The t-slope parameter Bp vs. W compared to predictions from the b-Sat and b-CGC 
models using the "boosted Gaussian" vector meson wave function. 
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Figure 22: The saturation scale Q 2 S = 2/r|, where r s is defined as the solution of found 
in the b-Sat and b-CGC models (left), and in the GBW and CGC models (right). 

that, with the exception of a' ¥ , the b-CGC model gives a considerably worse overall description 
of exclusive processes than the b-Sat model. 



5 Saturation and related topics 

A frequently asked question, whether or not the HERA data require saturation, is answered in 
the saturation models like GBW [35,36] or CGC [37] with a clear yes. In the impact parameter 
dependent models, such as the models discussed in this paper, the answer is more involved. In 
this section, we will therefore discuss the saturation effects in some detail. 

In the GBW model the effects of saturation are clearly seen, for example, in the change of 
rate of rise, A to t, of the total DIS cross section with Q 2 , see Fig. |S] (bottom). In this model 
the value of the observed parameter A tot is related to the value of the constant Agbw ~ 0.3 
modulated by the saturation effects. Since the variation of A tot with Q 2 is substantial, saturation 
has to be an important effect. Note that, in the GBW model, saturation is the only mechanism 
which can modulate the parameter Atot- 

The saturation effects are best quantified by the value of the saturation scale Q 2 S = 2/r|, 
where the saturation radius r s is defined as the solution of ()54|). In Fig. 1221 we show the 
saturation scale for the impact parameter dependent (left) and independent (right) models. 
Fig. 1221 (right) shows that the saturation scale in the GBW model is significantly higher than 
in the CGC model. (The GBW and CGC fits shown here are described in more detail in 
Sect. 15.21 ) This is understandable since in the CGC model the variation of the A tot parameter 
is partly due to evolution in addition to the saturation effects. However, even in the CGC 
model the saturation effects are fairly strong, as discussed by Iancu, Itakura and Munier [37] 
and by Forshaw and Shaw [38]. Fig. 1221 (left) shows that in the b-Sat and b-CGC models the 



35 



s-s 0.35 

© 0.3 

O 

w 0.25 

J3 

X5 0.2 

\ 

£ 0.15 

b 

ft 

b 0.05 

\ 

«- 

01 2345678 

b (GeV -1 ) 

Figure 23: The 6-dependence of the total cross section, a^, for Q 2 = 0.4, 4 and 40 GeV 2 with 
x = 10 -4 , 10 -3 and 10 -2 respectively, using a Gaussian T(b) of width = 4 GeV -2 . 

saturation scale is strongly dependent on the impact parameter b. In the centre of the proton 
(b ~ 0), the b-Sat and b-CGC models have a similar saturation scale, comparable to the value 
in the GBW model. As b increases the value of the saturation scale drops quickly in both 
models. This is again understandable since, in the b-Sat model with a Gaussian proton shape, 
at larger values of b the gluon density is diluted by the factor T(b) and so the smaller gluon 
density leads to smaller saturation scales. In this model, the variation of A to t with Q 2 is mostly 
due to evolution effects, since the gluon density at the initial scale /j/q is characterised by a low 
value of the parameter A 5 ~ 0. The observed large values of A to t can only be generated by 
evolution, as discussed in detail by KT [1]. 

In Fig.|2niwe show the 6-dependence of the total cross section to give a feeling for the relative 
contributions from the different impact parameters. The median value of this distribution is 
around b = 2.6 GeV -1 , that is, the majority of the cross section is determined by the dilute 
gluon region, where the saturation scale is small. 

To summarise, in the impact parameter dependent dipole models, evolution plays a greater 
role than saturation on average. However, in the centre of the proton (b m 0), the saturation 
effects are large in both the b-Sat and b-CGC models. In the centre of the proton the saturation 
scale is comparable to the saturation scale found in the original GBW model. 

The GBW model is theoretically very attractive since all observables in this model are a 
function of only one variable, r 2 Q 2 (x), where Q 2 (x) = (x /x) Agbw . This leads to so-called 
geometric scaling in which a 1 p is only a function of r = Q 2 /Q 2 {x), which is confirmed to 
some accuracy by data [51]. 8 A similar scaling has recently been observed for (^-integrated) 

8 Note, however, that the inclusion of the charm quark contribution violates geometric scaling to a certain 
extent. 
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diffractive DIS data [52]. The notion of geometric scaling is essential for development of the 
theoretical approach to saturation. Indeed, geometric scaling seems to be a universal feature 
of a wide class of evolution equations with saturation effects, irrespective of the form of the 
non-linear term [53-55]. In the b-Sat model, approximate geometric scaling is also present, as it 
is imposed by the fit to the data. This scaling, however, is not an intrinsic feature of the b-Sat 
model because of the greater importance of DGLAP evolution compared to saturation effects, 
and also because of the additional scale introduced by the impact parameter dependence. 

The theoretical understanding of saturation phenomena follows from evolution equations 
obtained within perturbative QCD. It is therefore interesting to ask the question whether the 
saturation effects determined in the models from fits to HERA data belong in the perturbative 
or non-perturbative domain. As shown in Fig. |22 the saturation scale determined in the 
proton centre in the b-Sat model is around 0.5 GeV 2 at x ~ 10 -3 . This number lies in- 
between the value of Aq CD = 0.04 GeV 2 , being clearly non-perturbative, and the value of 
around 1 GeV 2 , considered to be perturbative. Therefore it is not obvious to what extent 
the saturation dynamics are driven by the perturbative effects. The models discussed here 
are, however, by construction perturbative; the renormalisation and factorisation scale /i 2 = 
4/r 2 + //q, used to evaluate the strong coupling and the gluon density, is bounded from below 
by $ ~ 1 GeV 2 and is around 2 GeV 2 if Q 2 S = 2/r| ~ 0.5 GeV 2 . Moreover, in the centre of 
the proton, the value of the saturation exponent 

_ dln(Ql) 

Xs = dHl/xj (58) 

varies between As = 0.19 at x — 10 -2 and Xs = 0.27 at x — 10 -4 , as shown in Fig. EU 
Therefore, the values of this exponent are greater than the value of Xs — 0.08 expected for a 
'soft' process, and are close to the expectations from theoretical studies of perturbative non- 
linear evolution equations; see, for example, Refs. [56-58]. This indicates that the saturation 
phenomena studied in the b-Sat model is outside of the non-perturbative region. 



5.1 Step T(b) 

We also performed an alternative fit to F 2 data using the b-Sat model with the step function 
T{b) given by (}4*H|) . with the parameter bs = 4 GeV -1 ; see the last line of Table El Recall that 
this form of T(b) is implicitly used in all b- independent parameterisations of the dipole cross 
section. The fit was of similar quality and gave a slightly larger gluon distribution compared to 
the corresponding fit with a Gaussian T(b), see Fig. El indicating a slight shift in the balance 
between evolution and saturation effects. Note from ()49|) that a step T(b) with bs = 4 GeV -1 
corresponds to (b 2 ) = 8 GeV -2 , the same value as for the Gaussian T(b) with Bq = 4 GeV -2 
from (gBJ). 

For small \t\, the results with a step T(b) are close to those with a Gaussian T(b), and so 
the total cross sections for exclusive processes are also similar. However, for larger values of 
\t\, the step T(b) gives a dip in the t-distributions, which is not observed in the data, as seen 
in Fig. ESJfor J/tp production. Here, we have used the "boosted Gaussian" vector meson wave 
functions in both cases. The reason for the dip at large \t\ can be explained by noticing that 
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Figure 24: The saturation exponent \$ as a function of x and impact parameter b. 



Model 


Q 2 /GeV 2 


m U}d)S /GeV 


m c /GeV 


cr /mb 


xo/10" 4 


A 


X 2 /d.o.f. 


GBW 


[0.25,45] 


0.14 




20.1 


5.16 


0.289 


216.5/130 = 


1.67 


GBW 


[0.25,45] 


0.14 


1.4 


23.9 


1.11 


0.287 


204.9/130 = 


1.58 


GBW 


[0.25,650] 


0.14 


1.4 


22.5 


1.69 


0.317 


414.4/160 = 


2.59 


CGC 


[0.25,45] 


0.14 




25.8 


0.263 


0.252 


117.2/130 = 


0.90 


CGC 


[0.25,45] 


0.14 


1.4 


35.7 


0.00270 


0.177 


116.8/130 = 


0.90 


CGC 


[0.25,650] 


0.14 


1.4 


34.5 


0.00485 


0.188 


173.7/160 = 


1.09 



Table 5: Parameters of the GBW (J35|) and CGC p7jl models determined from fits to F 2 
data [22,23]. 



the two-dimensional Fourier transform of the step function (148)1 gives the Bessel function of the 
first kind, 2Ji(6sA)/(6gA), which oscillates through zero, whereas the two-dimensional Fourier 
transform of a Gaussian is simply another Gaussian. Although there is some uncertainty in the 
measured cross section at large |i| due to the treatment of proton dissociation, the uncertainty 
is not expected to account for the large discrepancy between the predictions with the step T(b) 
and the data, and so the step T(b) must be ruled out as a model for the proton shape. 

5.2 The GBW and CGC models without impact parameter depen- 
dence 

For completeness we give here the results of the fits using the impact parameter independent 
GBW flHHJ) and CGC dipole models. We first make fits to ZEUS F 2 data [22, 23] with 
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_i i i I i i i I i i i I i i i_jJ i i i I i i / 

0.2 0.4 0.6 0.8 1 1.2 



|t | (GeV 2 ) 

Figure 25: Differential J/ip meson cross section da/dt vs. \t\ compared to predictions from the 
b-Sat model using the "boosted Gaussian" vector meson wave function with either a Gaussian 
or a step T(b). 
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Figure 26: The saturation scale Q 2 S = 2/r|, where r s is defined as the solution of found 
in the GBW and CGC models with and without charm quarks. The presence of charm quarks 
dramatically lowers the saturation scale, especially for the CGC model. 

x < 0.01 and Q 2 e [0.25,45] GeV 2 using the CGC model (jSZj) with Af = 0.7 (fixed), first 
without any charm quark contribution as in the original paper [37], then including the charm 
contribution. We also show the effect of including the data with higher Q 2 > 45 GeV 2 . For 
comparison, we perform similar fits using the original GBW model (|35|). We take x = x B for 
light quarks and x = x B (l + 4m 2 /Q 2 ) for charm quarks. The light quark masses are taken to 
be m u> d, s = 0.14 GeV, with the charm quark mass m c = 1.4 GeV. The results of these fits are 
shown in Table El 

We note that the description of the data by the CGC model is sizably better than by 
the GBW model. This is presumably due to the lack of evolution effects in the GBW model 
and can be seen from the fact that the worsening of the x 2 value when the data points with 
Q 2 > 45 GeV 2 are included is more prominent for the GBW model than the CGC model; see 
the right-hand column of Table 

Notice also that the saturation scale in the CGC fit is dramatically lowered with the intro- 
duction of charm quarks, as shown in Fig. |2B The fact that saturation effects are very sensitive 
to the presence of the charm contribution was first noticed in the original GBW paper [35] and 
also in the KT [1] impact parameter dependent analysis. In particular, Thorne [59] has em- 
phasised the importance of the charm contribution, which has been omitted in some analyses 
of the saturation scale at HERA. 

In Fig. EH we show the dipole cross section a q q(x,r) at fixed x = 10~ 4 , integrated over the 
impact parameter b, obtained in the fits using the b-Sat, b-CGC, GBW and CGC models with 
m u,d,s = 0.14 GeV and m c = 1.4 GeV. At smaller values of r the b-Sat model has a slightly 
larger dipole cross section than the other models due to the presence of DGLAP evolution. At 
larger r the GBW and CGC models tend to a constant value of oo, while the b-Sat and b-CGC 
models continue to increase with increasing r; see also the discussion in Ref. [1]. However, as 
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r (GeV" 1 ) 

Figure 27: The dipole cross section a q g(x,r) at fixed x = 10~ 4 , integrated over the impact 
parameter b, obtained in the b-Sat, b-CGC, GBW and CGC models. 

discussed in Sect. 12. 1. lj. the contribution to the total cross section from large dipole sizes is 
generally suppressed by the photon wave functions, as is clearly seen in Fig. EHl 

6 Summary and outlook 

We have presented an analysis of exclusive diffractive vector meson and DVCS data measured 
at HERA within an impact parameter dependent saturated dipole ("b-Sat") model. Various 
cross sections measured as a function of Q 2 , W and t can be described by a model with a 
minimal number of free parameters, namely the parameters /Zq, A g and \ g of the initial gluon 
distribution, xg(x,pl) = A g x~ Xg (1 — x) 5 6 , and the proton width B G . The wave functions of 
the virtual photon are known from QED, while the vector meson wave functions are assumed 
to have a Gaussian shape. The variable which fluctuates in the Gaussian is, of course, not 
known precisely. However, we have shown that the observed distributions are fairly insensitive 
to the particular assumptions, with possible exception of the (JlI&t ratio for the p meson. A 
more precise measurement of this distribution and of the spin density matrix elements would 
allow better constraints to be made on the form of the vector meson wave functions. 

An important finding of this investigation is that, although the vector meson wave functions 
are not fully known, one obtains a good description of the measured data. The model param- 
eters, which were fixed by the fit to the total inclusive DIS cross section and the vector meson 
t-distributions, describe the measured Q 2 and W dependence of vector meson production and 
DVCS very well, together with the absolute normalisation. The measured DVCS t-distribution 
agrees with the model expectation within the measurement error. We expect that the high 
luminosity achieved by HERA will allow the t-distributions of vector mesons and DVCS to be 
measured more precisely. They provide important information about the proton size and the 
transverse dynamics of the evolution process. 
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Figure 28: The distribution of dipole sizes r contributing to the total inclusive DIS cross section 
in the b-Sat model for various virtualities, Q 2 , of the photon. The median values are indicated 
by vertical arrows. 
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The b-Sat model, which gives the best description of data, uses the Glauber-Mueller dipole 
cross section (J41j) with DGLAP evolution of the gluon density. Although the overall description 
of exclusive processes is very good, this approach has some limitations, seen most clearly in the 
lack of W dependence of Br, in J ftp photoproduction, Fig. [TBI Although this is a delicate effect, 
the measurement precision is sufficient to show that there is a coupling between the transverse 
and longitudinal evolution variables, that is, a' P ^ 0. We therefore introduced impact parameter 
dependence into the CGC model, the "b-CGC" model, which leads to a considerably poorer 
fit to F2 than the b-Sat model and a worse overall description of exclusive processes, but a 
better description of the a' P parameter. The saturation scale Q 2 S evaluated in this investigation 
does not depend sizably on the adopted evolution scheme and is consistent with the results of 
Ref. [1]. 

An important finding of this investigation is that the t-dependences of all three vector 
mesons and the DVCS process can be simultaneously described with one universal shape of the 
proton. The parameter characterising the size of the proton, Bq = 4 GeV -2 , determined in 
this investigation, corresponds to a root-mean-square impact parameter a/ (b 2 ), given by (JH5J), 
of 0.56 fm. This is rather smaller than the proton charge radius of 0.870 ± 0.008 fm [60]. 9 This 
leads to a rather surprising result that gluons are more concentrated in the centre of the proton 
than quarks. DVCS measurements planned at JLab should help clarify this somewhat puzzling 
picture (see, for example, [61]). 

The investigation presented here demonstrates that a wide class of high-energy scattering 
processes measured at HERA may be understood within a simple and unified framework. The 
key ingredient is the gluon density which is probed in the longitudinal and transverse directions. 
The success of the description indicates the universality of the emerging gluon distribution. 

Let us finish with a general remark that vector meson and DVCS processes may be used to 
probe the properties of nuclear matter in a new way. In measurements with polarised beams 
it is possible to achieve precision which would allow a holographic picture of protons and 
nuclei to be obtained [62-64]. Such a measurement could be performed at an ep collider with 
roughly a third of the HERA centre-of-mass energy, similar to the one described in the eRHIC 
proposal [65,66]. 
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A Connection to the KT paper 

In the preceding analysis [1] of J/ ip photoproduction in the impact parameter dependent dipole 
saturation model, Kowalski and Teaney (KT) used a somewhat different convention to define 

9 The proton charge radius was first measured by Hofstadter [48] to be 0.74 ± 0.24 fm. 
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the wave functions and to calculate the decay constants and the overlaps. 

KT [1] defined the overlap functions between the vector meson and the photon wave func- 
tions in the following way: 



(*^) T = e /e £^ {m}K (er)Mr,z) 



2-Kmf 
2tt 



(V* v V)l = e f e^-^2QK {er)z{\ - z )( p L {r, z„ 



[z 2 + (1 - zf] eK^drMr,*)} , (59) 

(60) 



where the scalar "Gaus-LC" wave functions 4>t,l{ t iZ) were defined as the Fourier transforms 
of factorised wave functions given in the momentum space by 



(j) TiL (k, z) = N T:L z(l - z) exp(-k 2 R 2 TL /2) 



(61) 



leading to 



d 2 k 



(2tt) 2 
N t ,lz(1 - z 



exp (ifc • r) (f>T,L(k, z 
d 2 k 



(2tt 



exp (ifc • r) exp(-A; 2 i?| L /2) 



2nR 2 TL 



z(l — z) exp 



2R\ L 



(62) 



In that representation the normalisation conditions were given by 



d 2 k r l dz 



(2vr) 2 



- { ^ 2 + (l-^) 2 ]^ + l!> MKz) 



m. 



d 2 k r x dz 



4> L {k,z) 



(63) 
(64) 



and the decay constants read, 



f v , T = e f V2N c —J—J o _ g) 



dz 



[z 2 + (1 - z) 2 ] " + 1 Mk 
nit 



(65) 
(66) 



It is straightforward to observe that the KT formulae (|59I6U|) . (|63I64|) and (|65I66|) may be 
obtained from the formulae of the present paper (|21|22j) . (|24I25|) and (|26I27|) if 5 = and the 

previously used wave functions 4>t and <pL are expressed in terms of the wave functions (f>T 
and <pL written in the conventions of this paper: 



ir(r, z) = - rrif <p T {r, z), 
h(r,z) = s/2N c M v 4> L {r,z). 



(67) 
(68) 
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Meson 


My/GeV 


fv 


m f /GeV 


Nr/GeV- 1 


R 2 T /GeV- 2 


N L /Ge\- 1 


R 2 L /GeV- 2 




3.097 


0.274 


1.4 


111 


6.5 


119 


3.0 




1.019 


0.076 


0.14 


164 


16.0 


214 


9.7 


p 


0.776 


0.156 


0.14 


211 


21.9 


222 


10.4 



Table 6: Parameters of the "Gaus-LC" vector meson wave functions. These are identical to 
those in Table [TJ but using N T>L instead of N TjL ; see (J6T?|) and (f7U|) . 



Note the modification of the ^-dependent part of (j)T(r,z). Of course, the radius parameters 
Rt,l are the same in both conventions. The normalisation factors are, however, transformed 
according to 

N T = \pm c m f 2tiR 2 t N t , (69) 
N L = jW c My 2nR 2 L N l . (70) 

The parameters of the "Gaus-LC" wave functions in its initial formulation are given in Table |H1 
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